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CN ' Abstract 

o : 

(30 \ We investigate the Q 2 evolution of parton distributions at small x values, recently 

Q\ \ obtained in the case of soft initial conditions. The results are in excellent agreement 

^ ■ with deep inelastic scattering experimental data from HERA. 

Oh! 

The measurements of the deep-inelastic scattering structure function F 2 in HERA [[TJ have 
permitted the access to a very interesting kinematical range for testing the theoretical ideas 
on the behavior of quarks and gluons carrying a very low fraction of momentum of the proton, 
the so-called small x region. In this limit one expects that non-perturbative effects may give 
essential contributions. However, the resonable agreement between HERA data and the NLO 
approximation of perturbative QCD that has been observed for Q 2 > lGeV 2 (see the recent 
review in J2|) indicates that perturbative QCD could describe the evolution of structure func- 
tions up to very low Q 2 values, traditionally explained by soft processes. It is of fundamental 
importance to find out the kinematical region where the well-established perturbative QCD 
formalism can be safely applied at small x. 

The standard program to study the small x behavior of quarks and gluons is carried out 
by comparison of data with the numerical solution of the DGLAP equations [] by fitting the 
parameters of the x profile of partons at some initial Ql and the QCD energy scale A (see, for 
example, ||). However, if one is interested in analyzing exclusively the small x region, there is 
the alternative of doing a simpler analysis by using some of the existing analytical solutions of 
DGLAP in the small x limit (see for review). This was done so in Ref. @]-|J where it was 
pointed out that the HERA small x data can be interpreted in terms of the so called doubled 
asymptotic scaling phenomenon related to the asymptotic behaviour of the DGLAP evolution 
discovered in M many years ago. Here we will illustrate results obtained recently in ||. 

1. Thus, our purpose is to demonstrate the small x asymptotic form of parton distributions 
in the framework of the DGLAP equation starting at some with the flat function: 

f a (Ql) = A a (a = q,g), (1) 



1 At small x there is another approach based on the BFKL equation, whose application is out of the scope 



where f a are the parton distributions multiplied by x and A a are unknown parameters that 
have to be determined from data. Through this work at small x we neglect the non-singlet 
quark component. 

In H an effective method to reproduce the x-dependence of parton distributions has been 
developed. It is based on a separation of the singular and regular parts of the exact solution 
for the moments of parton distributions. The method allows in simplest way to reproduce the 
leading order (LO) results [|J and to construct the ^-dependence of parton distributions at 
next-to-leading order (NLO): 
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f q (x,Q 2 ) + -fa(Q 2 )f g (x,Q 2 : 



where I v (cr) is modified Bessel function, s and p are given by s 
a{Ql) - a(Q 2 ), and 



ln(a(Ql)/a(Q 2 )), p 



a = 2y (d + s + D + p)lnx 
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Inx 2ln(l/x) 
are NLO generalizations of the corresponding Ball-Forte variables and 
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A: = A q - 20a(Q 2 )A+ 



(3) 



Wherever in this work we use the notation a = a s j (47r). 

The nonzero components of the singular and regular parts of the terms d± and D± have the 
form: 
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where /3 is the first coefficient of QCD /5-function, ( n are Euler (^-functions and / is the number 
of active quarks. 



Our NLO results coincide with the corresponding of Ball and Forte in Ref. || if one 
neglects the "— " component, expands our NLO singular terms (p) k Ik+\{cr) in the vicinity 
of the point a = o^o and ignores the NLO regular terms (i.e. put exp( — D + (l)p) = 
exp(—D_(l)p) = 1 and cancel the terms proportional to oc(Qq) into the normalization 
factors A^ and A~). We think, however, that this expansion is not so correct because it 
generates NLO corrections of the order of the LO terms. 

The behaviour of eqs. fl2|) can mimic a power law shape over a limited region of x, Q 2 : 

f a (x, Q 2 ) ~ x - x ° !{x > q2) and F 2 (x,Q 2 ) ~ x " a m (*>0 a ) 

The quark and gluon effective slopes \ e Jf are reduced by the NLO terms that leads to 
the decreasing of the gluon distribution at small x. For the quark case it is not the 
case, because the normalization factor of the "+" component produces an additional 
contribution undampening as ~ (Inx) . 

The gluon effective slope Xy* is larger than the quark slope Ag", which is in excellent 

agreement with a recent MRS and GRV analyses M. 

Indeed, the effective slopes have the asymptotical values (at large Q 2 ): 

K ff ' as (x,Q 2 ) = pW\~p 
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where symbol ~ marks approximations obtained by expansions of modified Bessel func- 
tions I n (a). The slope X e p2 ,as {x,Q 2 ) lies between quark and gluon ones but closely to 
quark slope X't^ ,as (x, Q 2 ) (see also Fig. 2). 

• Both slopes X e Jf decrease with decreasing x. Ax dependence of the slope should not 
appear for a PD with a Regge type asymptotic (x~ x ) and precise measurement of the 
slope X e Jl may lead to the possibility to verify the type of small x asymptotics of parton 
distributions. 

2. With the help of the results obtained in the previous section we have analyzed F 2 HERA 
data at small x from the HI collaboration (first article in In order to keep the analysis 
as simple as possible we have fixed A|^[g"( n / = 4) = 250 MeV which is a reasonable value 
extracted from the traditional (higher x) experiments. The initial scale of the PD was also 
fixed into the fits to Q% — 1 GeV 2 , although later it was released to study the sensitivity of 
the fit to the variation of this parameter. The analyzed data region was restricted to x < 0.01 
to remain within the kinematical range where our results are accurate. Finally, the number of 
active flavors was fixed to /=4. 

Fig. 1 shows F 2 calculated from the fit with Q 2 > 1 GeV 2 in comparison with HI data. 



closer to the data than the LO curve (dashed line). The lack of agreement between data and 
lines observed at the lowest x and Q 2 bins suggests that the flat behavior should occur at Q 2 
lower than 1 GeV 2 . In order to study this point we have done the analysis considering Qq as a 
free parameter. Comparing the results of the fits (see 0) one can notice the better agreement 
with the experiment of the NLO curve at fitted = 0.55GeV 2 (solid curve) is apparent at 
the lowest kinematical bins. 

Finally with the help of Eqs. @ we have estimated the F 2 effective slope using the value 
of the parameters extracted from NLO fits to data. For HI data we found 0.05 < A|4 < 
0.30 - 0.37. The lower (upper) limits corresponds Q to Q 2 = 1.5 GeV 2 (Q 2 = 400 GeV 2 ). The 
dispersion in some of the limits is due to the x dependence. Fig. 2 shows that the three types 
of asymptotical slopes have similar values, which are in very good agreement with HI data 
(presented also in Fig. 2). The NLO values of \p 2 as he between the quark and the gluon 
ones but closer to the quark slope \ e J^ as . These results are in excellent agreement with those 
obtained by others (see the review @ and references therein). 

3. Resume. We have shown that the results developed recently in || have quite simple 
form and reproduce many properties of parton distributions at small x, that have been known 
from global fits. 

We found very good agreement between our approach based on QCD at NLO approxima- 
tion and HERA data, as it has been observed earlier with other approaches (see the review 
f2j). Thus, the nonperturbative contributions as shadowing effects, higher twist effects and 
others seems to be quite small or seems to be canceled between them and/or with ln(l/x) 
terms containing by higher orders of perturbative theory. To clear up the correct contributions 
of nonperturbative dynamics and higher orders containing strong ln(l/x) terms, it is necessary 
more precise data and further efforts in developing of theoretical approaches. 
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